Topology-induced confined superfluidity in inhomogeneous arrays 
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We report the first study of the zero-temperature phase diagram of the Bose-Hubbard model 
on topologically inhomogeneous arrays. We show that the usual Mott-insulator and superfluid 
domains, in the paradigmatic case of the comb lattice, are separated by regions where the superfluid 
behaviour of the bosonic system is confined along the comb backbone. The existence of such confined 
superfluidity, arising from topological inhomogeneity, is proved by different analytical and numerical 
techniques which we extend to the case of inhomogeneous arrays. We also discuss the relevance of 
our results to real system exhibiting macroscopic phase coherence, such as coupled Bose condensates 
and Josephson arrays. 

PACS numbers: 05.30.Jp, 73.43.Nq, 74.81. Fa, 03.75.Lm 



I. INTRODUCTION 

The Bose-Hubbard (BH) Hamiltonian, describing 
bosons hopping across the sites of a discrete structure 
and originally introduced to model liquid He in con- 
fined geometrieSji proves successful in capturing the es- 
sential physics of a wide range of condensed-matter sys- 
tems. The best known examples are no doubt provided 
by Josephson-junction arrays (JJA)2^ and Bose-Einstein 
condensate (BEC) arrays^, which are the subject of a 
huge amount of both theoretical and experimental on- 
going investigations. The hallmark of such class of sys- 
tems is the presence of a superfluid phase as opposed to a 
(Mott) insulator phase. 6 . The theoretical studies hitherto 
carried out on such phase transition mostly focus on ho- 
mogeneous ambient lattices, and provide well-established 
numerical and analytical techniques. Homogeneous lat- 
tices are also the basis of the current experimental real- 
izations of systems belonging to the BH class. This is 
at least partly due to present technical constraints. For 
instance the optical techniques used to fragment BECs 
yield quite naturally homogeneous arrays. 

However, the striking progress in experimental tech- 
niques suggests the realization of inhomogeneous net- 
works to be at hand. Actually, JJAs can be engineered in 
nontrivial geometries with the only possible constraint of 
planarity. In this respect, interesting geometry-driven ef- 
fects are proposed in Refs. 0H, while the physics of a frac- 
tal JJA is experimentally studied in Ref. 0- As to BEC 
arrays, two very promising approaches for realizing inho- 
mogeneous topologies are provided by holographic optical 
trapsiSiii and magnetic microtrapsiSiiii^. In the latter 
case, ongoing efforts are aimed at reducing the spacing 
between individual microtraps, currently bound above a 
few yitm, in order to couple the condensates therein con- 
fined. 

The deep influence of topological inhomogeneities on 
thermodynamic properties of discrete boson systems is 
evidenced by the occurrence of unexpected features even 



in the absence of boson interaction. Indeed, a finite- 
temperature Bose-Einstein condensation can take place 
despite the low di mensio nality of the system. This is il- 
lustrated in Refs. Il5ll6l in the case of the square comb 
lattice, namely an array of linear chains (fingers) joined 
along a transverse direction (backbone) such as the one 
shown in the inset of Fig. ^ More precisely, inhomogene- 
ity induces a hidden band in the single-particle energy 
spectrum which is ultimately responsible for condensa- 
tion. 

On the other hand, the rich zero-temperature phase 
diagram of the BH model ensues from the competition 
between the on-site repulsion and the kinetic energy of 
the boson gas. In the light of this, a natural question 
arises as to the influence of topology on the physics of 
interacting bosons. In this respect we mention that the 
effect of the inhomogeneity arising from the superposition 
of a local on-site potential on an otherwise regular lattice 
has been recently addressed. In particular, the existence 
of local Mott domains induced by a parabolic confining 
potential has been eviden ced for BEC arrays in Refs. hj, 
while Refs. I22l2.3l2 ll'25l analyze the phase 
diagram on superlattices. 

Here, we consider inhomogeneities of purely topologi- 
cal (i.e. kinetic rather than potential) origin, focusing on 
the emblematic case of comb lattices, where the larger 
connectivity of the backbone is expected to act as a cat- 
alyst for superfluidity. Interestingly, the competition be- 
tween kinetic and boson-interaction energy causes the 
occurrence of an intermediate domain in the BH phase 
diagram. The usual Mott-insulator lobes are separated 
from the superfluid domain by a phase characterized by 
the localization of superfluidity in a narrow region sur- 
rounding the comb backbone, the rest of the structure 
exhibiting an unexpected insulator-like behavior. More 
precisely, we show that the local compressibility^ features 
an exponential decrease with increasing distance from 
the backbone. Note indeed that the topological inhomo- 
geneity of the structure requires a description in terms 
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of site-dependent quantities. These results, which — to 
the best of our knowledge — are the first concerning the 
influence of topology on the BH phase diagram, required 
the generalization of different numerical and analytical 
techniques^. The presence of confined superfluidity is 
first evidenced within a mean-field approach, and further 
confirmed by both a third order analytical strong cou- 
pling perturbative expansion (SCPE) and (population) 
quantum Monte-Carlo (QMC) simulations. 



II. THE BOSE-HUBBARD MODEL ON A 
GENERIC STRUCTURE 

The BH Hamiltonian, describing locally interacting 
bosons on a generic discrete structure consisting of M 
sites, is 
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where the operator ol 
at site j and n,j = at a 



+ (aj) creates (annihilates) a boson 
j counts the bosons sitting at site 



j. As to the parameters, U > accounts for the (on- 
site) repulsion among bosons, p, is the chemical potential 
and T is the hopping amplitude between adjacent sites, 
specified by the adjacency matrix A. This is a useful tool 
supplied by graph theoryjSI allowing an algebraic descrip- 
tion of the topology of a generic discrete structure. Its 
generic matrix element A^j is one if sites (h,j) are near- 
est neighbors and zero otherwise. In view of [N, H] = 0, 
where N = ^ ■ rij , Hamiltonian can be conveniently 
studied exploiting its block-diagonal structure. Since we 
are interested in the zero temperature phase diagram, 
hereafter (•) denotes the expectation value on the ground 
state of H. 

As we mentioned above, in the case of homogeneous 
topology the competition between on-site interaction and 
hopping gives rise to an interesting zero-temperature 
phase-diagram in the fj,/U-T/U plane, where two dif- 
ferent domains can be recognized. A superfluid phase, 
where the energy cost of adding or subtracting a bo- 
son to the system vanishes in the thermodynamic limit; 
An incompressible Mott-insulator phase, consisting of a 
series of adjacent lobes, where such operations cost a 
finite amount of energy and the filling / = N/M is 
pinned to an integer value. The Mott-supcrfluid tran- 
sition can be furthermore characterized by the compress- 
ibility k — dN/dfx, which is finite in the superfluid re- 
gion and vanishes within the Mott lobes. In the case 
of inhomogeneous systems, the possible effects of topol- 
ogy can be described in detail by the site-dependent local 
compressibility, 19 Kj — dpj/dp, where pj — (rij) is the 
local density of bosons. 

Owing to the enormous size of the Fock space, an 
exact solution of the model cannot be faced even for 
relatively small structures. However, the essential ele- 
ments of the Mott-superfluid transition can be captured 



resorting to different approximate schemes, such as mean- 
fieldj22i22i SCPEfSi the renormalization approach^ and 
QMC computations 



III. MEAN-FIELD APPROXIMATION 

The key point of the mean-field approach of Ref. |26 
consists in the approximation 
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allowing to recast Hamiltonian (JTJ as the sum of on-site 
Hamiltonians H « Ti. = . TLj. In the simple case of a d- 

— '3 J 

dimensional (translationally invariant) lattice, TLj is site 
independent and one is left with a single site problem 28 



H{a) = M[— n(n-l) 



fin - 2dT(a + a + )a + 2dTa 2 }, 



subject to the self consistency constraint a = (a), 34 
where the so-called superfluid parameter a can be consid- 
ered real without loss of generality. The phase diagram 
of the homogeneous case can be obtained numerically 28 
or even analyticallyiS&i 3 ^ 

For non-homogeneous structures 
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and the ground state of Ti has the form \ip) = 
&)j { a h}), where \j\ {ah}} is the ground state of Hj. 
Thus the problem is solved by finding the set of real quan- 
tities {ctf l }ff- 1 such that 



(j;{a h }\aj\j;{a h }) = 



(5) 



This can be easily done numerically by means of self- 
consistent iterative algorithmiii 2 ^ Despite the approxi- 
mation in Eq. J5| strongly suppresses spatial correlation, 
some topological information is retained in the above 
mean-field formulation owing to the presence of the ad- 
jacency matrix A in Hamiltonian In this case pj, kj 
and aj are in general site-dependent quantities. For comb 
lattices these quantities are constant along the backbone 
direction, owing to the symmetry of the system. 

Fig. ^ shows the numerically determined mean-field 
phase diagram for a 100 x 100 comb lattice. In the regions 
1/ (where / is the integer filling) Kj = for all j's and 
the total number of bosons is pinned at N — JM. The 
system is therefore an incompressible Mott insulator. In 
the regions II/, Kj is finite, yet it vanishes exponentially 
along the fingers. The same behaviour is observed for pj, 
which is exponentially close to / with increasing distance 
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FIG. 1: Inset: an example of square comb lattice featuring 
11 sites both on the backbone and on the ribs—; Main plot: 
Mean-field Phase diagram of a 100x100 comb lattice. Differ- 
ent shades of grey denote different phases. The dashed line 
at T/U — 0.055 signals the set of parameters considered in 

Fig.m 
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FIG. 2: Behaviour of the local density of bosons pj for sites 
j along a finger of the comb lattice (j = backbone). The 
figure refers to a fixed value of the hopping amplitude T/U = 
0.055 and a finite interval of p/U (dashed line in Fig. EJ. 
The white density profiles correspond to the borders of the 
different regions of the phase diagram. As we discussed in the 
text pj = 1 inside region Ii, pj — > / in regions II/, whereas 
Pj tends to a not necessarily integer number in region III. 



from the backbone. Hence in these regions the superfluid 
behaviour of the system is confined along the backbone 
direction alone. An extended superfluid behaviour is re- 
covered in region III, where the local density pj far from 
the backbone is not necessarily an integer quantity, and 
Kj is nowhere vanishing. Such behaviour of the local 
density of bosons and compressibility are summarized in 
Fig. [21 and in the upper panel of Fig. |3 respectively. 

We mention that it is possible to evaluate the exact 
analytical form of the boundaries of the 1/ regions as pro- 
vided by the mean- field approach described by Eqs. Q- 
|[5J| . This can be for instance accomplished making use of 
the finite-temperature method reported in Ref. l25l and 
subsequently letting the temperature go to zero. The 
function of T/U describing the If boundary for a generic 
structure characterized by the adjacency matrix A is 
obtained by rescaling the corresponding function for a 
homogeneous <i-dimensional lattic o 25 i 35 by a factor A, 
where A is the maximal eigenvalue of A. 

In general, topological inhomogcncitics make the study 
of critical behaviours a rather difficult task. However, 
the above results suggest some considerations in this re- 
spect. On regular lattices the correlation length diverges 
in any direction at a critical point. Conversely, on the 
comb lattice, the correlation length is expected to di- 
verge only along a specific direction, depending on the 
critical border under concern. To wit, the correlation 
length between sites of the same finger is the only diver- 
gent quantity at the II / III transition, while it is finite at 
the border between regions If and II/ (where the diver- 
gent quantity is the correlation length between sites of 
different fingers). We also mention that preliminary re- 



sults based on the mean- field approach of Ref. l25l indicate 
that the above picture is robust at small finite temper- 
atures, and hence in principle accessible to experiments. 
In this respect we note that the three different phases in 
Fig. can be probed as in Ref. la, provided that the sys- 
tem is realized in terms of coupled BECs, possibly using 
holographic traps— Indeed, after the trapping poten- 
tial is released, the expanding atomic clouds should pro- 
duce either a one-dimensional 37 or a two-dimensional 3 ^ 
interference pattern depending on whether superfluidity 
is confined along the backbone or extended on the entire 
comb. 



IV. BEYOND MEAN-FIELD 

A step beyond the mean-field approximation consists 
in the strong coupling perturbative expansion. Indeed, 
time-independent perturbation theory in the hopping pa- 
rameter allows to obtain an analytical approximation 
of the Mott lobes, which, in the case of the linear 
chain, proves to be quite satisfactory already at the third 
orderi 3 ^ This approach, introduced in Ref. |3JJ for homo- 
geneous bipartite structures, is extended to any structure 
in Ref. |2(J Quite interestingly, it turns out that topolog- 
ical inhomogeneity gives rise to a third-order correction 
featuring an unusual dependence on the adjacency matrix 
describing the topology of the structure. Indeed, unlike 
the previously reported perturbative terms, depending 
only on the maximal cigenpair of A, the "topological cor- 
rection" depends on the entire spectrum of the adjacency 
matrix. The solid line in the inset of Fig.0|is the border of 
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FIG. 3: Local compressibility Kj for sites j along one finger of 
the comb lattice (j = is on the backbone). Dashed, solid and 
dotted lines refer to region Ho, Hi and III, respectively. Note 
the exponential decrease with increasing distance from the 
backbone characterizing the first two curves. Upper panel: 
Mean-field result. Lower panel: QMC data for a 12x12 
lattice. 



the Mott lobe Ii for a comb lattice as provided by the an- 
alytical third-order strong-coupling perturbative expan- 
sion reported in Ref.|2^. The above described exponential 
localization of superfluidity characterizing phases II/ is 
also captured by SCPE even at order zero. Indeed it can 
be easily shown that pj = f + C\vj\ 2 and kj = K\vj\ 2 , 
where C and K are normalizing constants and Vj is the 
j-th component of the maximal eigenvector of the adja- 
cency matrix, depending only on the distance dj from 
the backboneji^*i& Vj = exp[— d/asinh(l)]. First order 
SCPE confirms this behaviour, though with a different 
decay rate, as well as the above considerations about the 
correlation functions^ 

The data points in the inset of Fig. 0] satisfactorily 
agreeing with the perturbative curve, have been obtained 
with a population QMC approach^ adopting the resam- 
pling procedure described in Ref. ^3 Based on a gen- 
eralization of the power method for finding the maximal 
eigenpair of a matrix, this technique essentially amounts 
to a stocastic evaluation of the ground state of Hamilto- 
nian Q and therefore allows to study the zero tempera- 
ture phase diagram of the BH model. Note that both in 
the SCPE and in the QMC approach p is evaluated as 
the difference between the ground state energies of sys- 
tems whose total number of bosons differs by onem&2£ In 
particular, the border of phases 1/ are obtained consid- 
ering the energy cost of adding or subtracting one boson 
from the integer filling situation N = Alf. 

The QMC approach also confirms the existence of the 
confined superfluid phases II / . Fig. 0] clearly shows the 
transition between phases Ho and III. Indeed in the for- 
mer region the local density pj at a site far from the back- 
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FIG. 4: Inset: Mott lobe L for a comb lattice according 
to third order SCPE (solid line) and population QMC (er- 
rorbars). The latter refer to a 12x12 lattice. Main plot: 
QMC results for the on-site density of bosons on a comb lat- 
tice at T/U = 0.01. Upward triangles: average result for a 
site on the backbone. Downward triangles: average result for 
the farthest site from the backbone. Open and filled symbols 
refer to a 12x12 and a 16x16 comb lattice, respectively. The 
larger error on the QMC data (abscissae) is smaller than the 
symbol size. 



bone (downward triangles) is very close to / = 0, and it 
features a sudden increase only after entering phase III. 
Conversely, the local density on the backbone (upward 
triangles) is fractional also in region Ho, thus confirming 
that the superfluid is localized there. 

A further confirmation of confined superfluidity is pro- 
vided by the local compressibility profiles appearing in 
the lower panel of Fig. |21 obtained by means of QMC 
simulations in the case of a 12x12 comb lattice. In- 
deed, in the extended superfluid region the local com- 
pressibility is everywhere significantly larger than zero, 
whereas the curves relevant to the regions Ho and Hi 
feature a sharp decrease with increasing distance from 
the backbone. Note that the local compressibility within 
the Mott lobes 1/ is everywhere zero, since, by definition, 
the ground state of the system can be changed only if the 
chemical potential p is varied of an amount sufficient to 
cross the lobe border. 



V. CONCLUSIONS 

In summary, we reported the first analysis of the in- 
fluence of topological inhomogeneity on the phase dia- 
gram of interacting bosons, considering the emblematic 
case of comb lattices. This supplies a basis and a con- 
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ceptual framework for a more general study aimed at a 
deeper understanding of the role of topology in quan- 
tum phase transitions. Furthermore we suggested a pos- 
sible experimental setup, based on BEC arrays trapped in 
holographic potentialsji2*ii where the intermediate phase 
occurring on a comb lattice could be observedi£i2L2£ 

The recently disclosed relation between critical be- 
haviour and system-state entanglement^! provides a fur- 
ther context where the influence of geometry might play 
a significant role. In particular, inhomogeneous arrays 
have been recently proposed as quantum-information- 
processing devices^ We point out that, owing to the for- 
mal mathematical analogy between Heisenberg and BH 
models (see e.g. Refs. I2I andElt) . the results herewith pre- 
sented have also relevant implications for quantum spin 
systems on inhomogeneous structures. As a concluding 



remark, we observe that a comb lattice can be obtained 
by joining ID structures or appropriately removing the 
exceeding links from a 2D regular array£ This makes the 
structures considered here not only interesting from the 
theoretical point of view, but also very promising for ac- 
tual realizations based on JJA technology. 3 9 44 
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